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Abstract 

We consider the Hamiltonian structure of reduced fluid models obtained from a kinetic description of colli- 
sionlcss plasmas by Vlasov-Maxwell equations. We investigate the possibility of finding Poisson subalgebras 
associated with fluid models starting from the Vlasov-Maxwell Poisson algebra. In this way we show that 
the only possible Poisson subalgebra involves the moments of zeroth and first order of the Vlasov distribu- 
tion, meaning the fluid density and the fluid velocity. We find that the bracket derived in [Phys. Rev. Lett. 
93, 175002 (2004)] which involves moments of order 2 is not a Poisson bracket since it does not satisfy the 
Jacobi identity. 



1. Introduction 

In the kinetic framework, the dynamics of collisionless plasmas is provided by the Vlasov-Maxwell equa- 
tions which are dynamical equations for the phase space density /(x, v) of the charged particles (also called 
Vlasov density) and the electromagnetic fields E(x) and B(x) where x and v belong to M 3 . Here, we work 
with only one species of charged particles of unit mass and charge e for the sake of simplicity, but the 
generalization to several species is straightforward (see, e.g., Ref. |l|). The dynamical equations are 

/ = -v V/-e(E + vxB)-9 v /, (I) 

E = V x B - e J d 3 vfv, (2) 

B = -V x E, (3) 

where the dot designates the time derivative. Given nowadays computer capability, integrating numerically 
this kinetic model is too demanding for realistic laboratory plasmas. Moreover the dynamics is not easily 
analyzed in such a kinetic framework since it provides dynamical information at temporal and spatial 
scales which might be irrelevant. In addition, the plasma dynamics is more conveniently analyzed in the 
configuration space rather than in the phase space of the charged particles. As a consequence, there is a 
need for reduction in order to eliminate the irrelevant parts of these equations so as to obtain a much more 
tractable model. One way of doing this is to consider a fluid reduction which is obtained by considering the 
moments of the Vlasov distribution, e.g., the fluid density, the fluid momentum, the pressure tensor, etc... 

Starting from the dynamics of the zeroth order moment, a series of dynamical equations for the higher 
order moments is built. This is a priori an infinite set of differential equations. The reduction is obtained 
by truncating this set and closing the system of equations through a closure assumption, often based on 
a collisional argument or on a local thermodynamical equilibrium. According to the physical situation of 
interest, the choice of an order of truncation and of a closure assumption brings us a corresponding fluid 
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model. For instance, starting from the parent model, i.e., the Vlasov-Maxwell equations, many physically 
interesting fluid models have been derived (see, e.g., Refs. A particularly important property 

of the parent model is that it possesses a Hamiltonian structure, i.e., the equations of motion can be rewritten 
using a Hamiltonian functional H and a Poisson bracket {•, •} as F — {F, H}. Truncating and closing the set 
of equations might not conserve this property. This introduces dissipative or non-Hamiltonian terms in the 
equations of motion whether they have a physical origin (as described by phenomenological constants like 
diffusion constant, magnetic diffusivity or kinematic viscosity) or not. If these terms do not have a physical 
origin, they have been coined mutilations in Ref. Q. If one considers only the ideal part of the equations 
which is the part where the phenomenological constants characterizing the dissipative terms are set to zero, 
one should recover a Hamiltonian system as an inheritance of the Hamiltonian parent model [H, HlSEll 



There are two ways to proceed. The first one is to work with the equations of motion and check a 
posteriori that the resulting set of reduced equations possesses a Hamiltonian structure (by finding an 
appropriate conserved quantity and a Poisson bracket). This can be tedious as one needs to check that 
the bracket which has been guessed satisfies all the properties of a Poisson bracket, and in particular the 
Jacobi identity. The second method is to work directly on the Hamiltonian structure of the parent model 
by performing the reduction on the Hamiltonian and on the bracket. Of course the main difficulty resides 
in verifying that the reduced bracket still satisfies the Jacobi identity. This puts some restrictions on what 
is allowed to do on a Poisson bracket in the course of the reduction. The gain is significant since the 
preservation of the Hamiltonian structure is ensured, i.e., there is no fake dissipation or mutilation, and 
it allows to keep track of the conserved quantities throughout the derivation and apply all the techniques 
already available for Hamiltonian systems (like perturbation theory, energy-Casimir methods for equilibria, 
Lie transforms, etc.). 

In Sec. [H we consider a very practical and useful method for deriving reduced fluid models by using 
Poisson subalgebra. We apply it to the Vlasov-Maxwell Hamiltonian structure. This method provides the 
Hamiltonian structure of the usual fluid model composed by the continuity equation for the fluid density 
and the momentum equation for the fluid velocity. However we show in Sec. [3] that when the closure occurs 
at orders higher than one, the reduced structure is no longer a Poisson subalgebra associated with the 
parent structure. As a consequence it is not possible to obtain reduced fluid models containing higher order 
moments by just applying a Poisson subalgebra argument. 

In Refs. a fluid model for low temperature relativistic plasmas, called the warm fluid model, 

involving second order moments of the Vlasov distribution was proposed. It aims at investigating the 
interaction between a strong laser pulse and a low-density plasma. In Ref. (l2j a conserved quantity has 
been constructed based on the Vlasov-Maxwell Hamiltonian, and a bracket was proposed. It has been 
shown that this model conserves energy and entropy. However, based on the results explained in Sec. [2J 
the Hamiltonian property of this model has to be scrutinized since it involves higher order moments of the 
Vlasov distribution. In Sec. [4] we exhibit a counter example for the Jacobi identity for the bracket proposed 
in Ref. [12j . Therefore, this warm fluid model is not a Hamiltonian system. 

2. Reduction to the usual fluid model 



The starting point is the Hamiltonian structure of the Vlasov-Maxwell equations 0, El El, EH- The 
Hamiltonian functional is the total energy : 



ff(/,E,B) = J d 3 vd 3 xf^- + J d 



, E 2 + B 2 
x . 



The Poisson bracket acts on the Poisson algebra of observables, that is the set of functionals of the field 
variables /(x,v), E(x) and B(x) : 

{F, G} v = J d 3 xd 3 v f[Ff,G f ] +ej d 3 xd 3 v f (d v F f ■ G E - F E ■ d v G f ) 

+ |A (F E -VxG B -f B -Vx G E ) , (4) 



where indicates the functional derivative with respect to the field variable ip and the bracket [■, •] is given 
by 

[/, g] = V/ ■ d v g - d v f ■ V.g + eB ■ d v f x d v g, 

and d v designates the partial derivatives with respect to v. The dynamics of a functional F of the Poisson 
algebra is given by F = {F, H}. In particular we recover Eqs. CO)-© if we choose F — /, F — E or F = B. 
We notice that we have considered here a non-relativistic plasma. However the discussion which follows is 
unchanged in the case of a relativistic plasma since the changes only occur in the Hamiltonian, not in the 
Poisson bracket [l]} . 

Fluid models rely on the idea of replacing the Vlasov density / by the series of its moments [19[ 



P^-'" = J d 3 vfv il ---v in , (5) 

for n G N. The chain rule yields 

oo 

Ff = }^F Pn i li2 ...i n u^Uij • • -v in , 

where we have used Einstein's convention of implicit summation over repeated indices ik- The expression of 
the Poisson bracket in these new variables involves the derivatives with respect to v given by 

oo n 
n=0 fe=l 

We use the following symmetrization of the tensor P n 

n-l 



FpiH- 'in) - 

k=0 



1 " ^ 

fc=0 

In such a way, the derivative with respect to v becomes 

oo 

d "i F f = ^2 n F p (si~-in) v n ■ •■v in _ 1 6\ n . 

n=0 

In these variables, the Hamiltonian becomes 

/Tyii f Tp2 _r_ T32 

d 3 x^- + J d 3 x — ±— , (6) 

and the expression of the bracket becomes 

{F,G} V = ^ f $ x mP n+m-l 9 kFp^G p (p k ) + ^2 I d ^ x nmP ni m -2—^- F p^ G P W) 

7n,n m,n 

+ Y. I ' d3x enP^_ 1 F p ^ ) G Ez + f d 3 x F E ■ V x G B - (F O G), (7) 

where Bij = £ijkB k with ey*. the Levi-Civita tensor, and (F -H- G) indicates that the terms obtained by 
inverting F and G in the summation have to be subtracted (in order to fulfill the antisymmetry property of 
the Poisson bracket). Here and in what follows, dk designates the partial derivative with respect to Xk] it 
acts on the next term, e.g. dkfg = (dkf)g- The greek indices a and f3 in Eq. ([7]) denote a set of indices so as 
to complete the summation. For instance, in P« , the indices a is a set of n — 1 indices a = . . . , i n -i) 
so thatpl m) =P ( l il - l "- ll) . 
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The commonly used fluid model corresponds to a truncation of the moments at order one, keeping as 
field variables only the fluid density p = Pq and the momentum density M = Pi (or equivalently the fluid 
velocity defined as M/p). The bracket ([7]) has the particular property that the set of all functionals of the 
reduced field variables (p, M, E, B) is invariant, in the sense that the Poisson bracket ([7]) of two functionals 
F(p, M,E,B) and G(p, M, E, B) is again a functional of (p, M, E, B). For this subset of functionals, the 
bracket reduces to 



{F, G}i 



d 3 x ((P d k F Po + PfdkFps)Gp* + eP Q (ByFptGps/2 + F P} G El ) + F E • V x G B 
(FoG), 



(8) 



which is indeed again a functional of (p, M, E, B). At this stage, the reduction of the dynamics has not 
yet been performed because the Hamiltonian H given by Eq. ([6]) depends on P2 so it does not belong 
to the subalgebra of functionals of (p, M, E, B). In order to perform the reduction, one needs to make 
an assumption on the Hamiltonian. For example, in the so-called cold plasma model [2J, the reduced 
Hamiltonian is H* = J d 3 x (M 2 /(2p) + (E 2 + B 2 )/2), i.e., the original kinetic part of the Hamiltonian 
/ d 3 x P"/2 has been replaced by J d 3 x M 2 /(2p). Practically, the choice is inspired by the physics of 
the system under consideration. In the presence of scalar pressure terms and in the absence of heat flux 
(adiabatic closure), the entropy s has to be included following Ref. To account for this, s has to be 
considered as an independent scalar field advected with the fluid. This is achieved by including the advection 
term given by 

d 3 x G S F M ■ Vs - (F o G), 



into the bracket. The addition of these terms preserves the Jacobi identity, since it only corresponds 
to a passive scalar advected by the Hamiltonian fluid. The last step is to include the pressure term in 
the Hamiltonian which may be any function pU{p, s). Accordingly, U(p, s) stands for the thermodynamical 
internal energy per unit mass, and the pressure is p 2 d p U. The model we get is the most common Hamiltonian 
fluid model (see Ref. (l6() with the Hamiltonian 



d 3 x ( — +pU(p lS ) + — 



B 1 



and the bracket 



{F,G} = 



d 3 2 



p Gm • VF P + M • (G M • V)F M + G s F M ■ V s + epB ■ F M x G M /2 

+epF M ■ G E + F E ■ V x G B ] - (F o G). 

The introduction of U in the Hamiltonian is a way of keeping some information that is enclosed in the 
pressure tensor. As a result, this model is almost a reduction at order 1 in the moment series, but it is 
related to a closure at order 2 as well. 



3. Fluid reduction at orders higher than 1 

As it was shown above, the Hamiltonian fluid reduction at order 1 is given by a Poisson subalgebra of 
the parent Hamiltonian structure. This results in a reduced bracket that is just a truncation of the initial 
bracket l[7|). However, in many cases such a fluid model with only the first two moments is not rich enough to 
account for all the physics of interest, and one has to include kinetic effects and retain higher order moments 
in the fluid model. In this section we consider the Hamiltonian derivation of such higher order models. It 
should be pointed out that including the second order moments would be particularly interesting since the 
Hamiltonian belongs to the reduced algebra, and the closure does not affect the Hamiltonian. 

The goal is to consider generalizing the results of the previous section to get fluid models of order higher 
than 1 by truncating the initial bracket. This method seems natural in the case of the fluid reduction, which 
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relies on the idea that all the physical information is contained in the set {P n }n<iV of the N first moments 
of /. The remaining moments will have to be expressed as functionals of them, Pi>jv — &i{{P n }m^N)- Here 
we consider the simplest closure assumption P n = for n > N. 

It should be noticed that it would actually be best suited to consider moments defined by quantities like 
n„ = J d 3 /(v — M/ p)® n I p. However it makes the discussion slightly more complicated, so for the sake 
of clarity, we have kept the moments as defined by Eq. ([5]) in this section. We consider reduced brackets 
defined from the moments Il„ in Sec. [H 

We look for the reduced bracket by removing the undesired moments from the initial Poisson bracket (|7J), 
i.e., we consider the sub-algebra of functionals F(f, E, B) = F({P n } n <N, E, B) and we truncate the Poisson 
bracket |7J) acting on two functionals of this sub-algebra by removing all terms proportional to P n for n > N . 
This bracket does not involve functional derivatives Fp n for n > N since these terms vanish when acting of 
an element of the subalgebra. We denote this truncated bracket {F, G}n- Since {F, G}n = {F, G}y for all 
functionals F and G of this sub-algebra by enforcing P n — for n > N, the reduced bracket is automatically 
algebraically closed and it retains from the initial bracket {•, -}y the bilinearity, the antisymmetry and the 
Leibniz rule. The only property one has to check to have a Hamiltonian structure is the Jacobi identity, 
since this property does not automatically transfer to truncated brackets. 

In order to be more specific, we inspect more closely the case N = 2. The reduced bracket is given by 

{F,G} 2 = {F,G} 1 + {F,G}' 2 , 

where the bracket {F, G}i is given by Eq. © and {F, G}' 2 is a bracket which only involves Fp 2 or Gp 2 . This 
bracket is given by 

{F, G}' 2 = J d 3 x (2(P(d k F Po + P; j d k F P} )G p u k) + l'?ih,F r G Pf + 2eP}F p w G Ej 

+2eB ij (p£F p »»G p ,+Pf l F p i ik) G p ui ) )) - (F o G). 

The bracket given by Eq. ([7]) between functionals of (p, M, P2,E, B) involves only one term that is not a 
functional of (p, M, P 2 , E, B), and which is proportional to P3, and it is given by 

{F,G}% =2 Jd 3 x P; jk diF p iiG p m -(F« G). 

The reduction involves P3 = 0, it is why this term has been dropped from the reduced bracket {•, - } 2 . 
However we can not conclude that the truncated bracket {•, -} 2 is Hamiltonian since the Jacobi identity has 
to be checked a posteriori. We show that in fact some terms not proportional to P3 are generated in the 
Jacobi identity by the contribution {•, - j^'. These contributions originate from the terms in the bracket ((7]) 
that are proportional to P n for n < 2 and which involve Fp 3 or Gp 3 . These terms are 

{F,G}° 3 = 3 Jd 3 x (Pl°d k F Po G p ^ k) + eB l3 P^ l F P} G p u kl) + eP^ F p « m G Eh ) - (F o G) 

If we restrict this additional contributions to the first term, i.e., in absence of magnetic and electric field, the 
non-trivial contribution in the Jacobi identity comes from the bracket {{F, G} 2 , which includes terms 
of the form 

-2 J d 3 x ptidiFptiGpwdkHpt, 

and as well two other terms, that have the same expression, but with circular permutation of (ijk). This 
contribution to the Jacobi identity of the bracket ([7J suggests a counter-example for the failure of the Jacobi 
identity for the truncated bracket {•, -} 2 . We select F as a functional of P 2 , G as a functional of P 2 and 
H as a functional of Pq. For instance, we choose F = J d 3 x (P 2 11 ) 2 /2, G = J d 3 x P 2 22 and H = Pq as a 
counterexample of the Jacobi identity. It leads to 

{{F, G} 2 , H} 2 + 0= -d\ ((P 2 n ) 2 ) - Ad, (P 2 12 <9 2 P 2 n ) . 
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The same reasoning can be performed at order TV. The number of families of possible counter-examples to 
the Jacobi identity increases as N increases. However a convenient one is inspired from the case N = 2 
and involves a functional of Pjv, a functional of P 2 and a functional of Pq. It follows that the truncated 
bracket {•, -}iv neglects terms proportional to Pjv+i which provides a non- vanishing contribution to the 
Jacobi identity for {•, -}n even in the subalgebra of functionals of {P n } n <N, E and B. 

Therefore, starting from N = 2, the truncated brackets {•, -}n do not satisfy the Jacobi identity. 

We notice that the reduced bracket {•, - }n contains most of the terms of the initial bracket {•, - }v when 
acting on functionals of the reduced variables {P n }n<N, E and B. It is then more efficient to study the few 
removed terms than the many kept ones. This suggests a simpler way of verifying the Jacobi identity by 
only considering the truncated terms. It is explained in the Appendix. 



4. Brackets expressed using pressure-like moments 



In Refs. [12. lla. 6], a warm fluid model for a collisionless low temperature relativistic plasma has been 
introduced. It gives dynamical equations for moments of the kinetic (Vlasov) distribution up to order two, 
i.e., it involves P (fluid density), IT = Px/Pq (fluid velocity) and 



It is argued in Refs. [12], |6[ that this model is Hamiltonian and the bracket is given by 

{F,G} 2 = {F,G} 1 + {F,G}' 2 , (9) 

where 

{F, G}i - J d 3 x (d k F Po G U k + ^(d k n[ - + eB M )F n uG nli + eF n *G Bfc + F E ■ V x Gr\ —(F O G). 



and 



{F,G}> 2 = J d 3 x (^^F nk G nt)+2 U^d k G <k) 

nr\ rs \ 
+— Z-(d k lj[ - dtUl + eB kl )F (kr) G n(ls) - (F O G). 

Pq 11 2 11 2 J 

This model is derived in a very similar way as in the previous section, except that, instead of going from 
/ to P n , the change of variables is from / to n„ where IIo = Pa, II 1 = Pi/Pq, and for n ^ 2, n„ = 
J d 3 v (v — ni)®"/Po- The reduction at order N corresponds to setting II n = for n > N. In particular, 
in the derivation of the model given in Ref. (l2l | , the contributions in the bracket which are proportional to 
II3 have been neglected given a specific assumption on the Vlasov distribution. However, as in the previous 
section, these neglected terms contribute to satisfying the Jacobi identity. Given the counter-example 



F = 


P^l 1 




G = 


jd* X 


p nf 


H = 




p nf 



the bracket @ fails to satisfy the Jacobi identity since 

{{F,G} 7 H}+ 0= 8P Ul 2 d 2 (^ 3 (p ni 3 )) +4d 2 ul 1 d 3 (p n 2 2 3 ) - (2 o 3), 

where (2 <-> 3) indicates the same two terms where n^ 2 has been exchanged with n3, 3 , and d 2 with 83 
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5. Conclusion 

In Ref. [l2j], it is mentioned that the reduction "is exact: for any functionals F[n, P, II] and G[rt,P,II] 
we have {F, G}m — {F, G}v- As a consequence we see that the moment bracket inherits the Jacobi identity 
(as well as all other properties) from the full bracket." Here the bracket {•, -}m refers to the reduced bracket 
{-, -}2 given by Eq. © and {■, -}y refers to the Vlasov-Maxwell bracket The first part of the sentence is 
correct. This is why the reduced bracket inherits all the properties of the Vlasov-Maxwell bracket that are 
linked to the values of the bracket: bilinearity, antisymmetry and Leibnitz rule. However the Jacobi identity 
involves not only the values of the bracket but also their gradients, and in general, 

{{F, G} M , H} m + {{F, G} v , H} v , 

even if {F, G}m = {F, G}v (after some assumption on the Vlasov distribution which is not preserved by 
the flow, or after neglecting higher order terms), so the reduced bracket can fail the Jacobi identity. 

Here we have shown that deriving reduced fluid models from kinetic models like Vlasov-Maxwell equations 
has to be handled with care as soon as moments of order two or higher are considered. As an illustration 
we have shown that the model introduced in Ref. [l2| involving moments up to order two does not satisfy 
the Jacobi identity, by exhibiting a counter-example. 
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Appendix: Verifying the Jacobi identity for truncated brackets 

First we consider a Poisson bracket written as 

{F,G}=F Zi F iZi G Zj . 

The Poisson matrix J satisfies the Jacobi identity rewritten as J2(ijk) § ZiZl d Zl I ZjZk = 0, where Ylajk) means 
circular permutation of the indices i,j,k. For the reduction we split the set of variables in two subsets 
z = (Z, £). The reduction is chosen as acting on the subalgebra of functions of Z. In order to do that, we 
decompose the Poisson matrix into 

j = j zz + sf z + ii z + jp 1 + j£ z] + j zz + ji z , 

where [£z] means antisymmetrization, i.e., JJ^ Z 1 = J^ z — JJ Z ^, the index in Jj indicates the i-th order in and 
Jr includes the remaining (higher order) terms. The reduced bracket is simply the truncation of Jf obtained 
by removing terms involving and taking the limit £ — > 0. As a result, 

JJ* = Jq 

This matrix J* is associated with a Hamiltonian system if and only if 

]T j z ' Ci a Ci jf^ fc =o. 

(ijk) 

In the example of Sec. [31 J 2 ** corresponds to {•, •}§, and J zz to {•, - j^. We notice that a Lie subalgebra 
corresponds to the case J zz = and JJ ZZ = 0. The condition on the validity of the reduction can be 
generalized to the case J zz = which is sufficient to verify the Jacobi identity. However this assumption is 
not often verified since brackets generally contain linear terms in the variables [19j . 
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